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Abstract. We propose an action-based f(R) modification of Einstein's gravity that 
\l admits a modified Schwarzschild-deSitter metric. In the weak field limit this amounts to 
C^l adding a small logarithmic correction to the newtonian potential. A test star moving in 

o '. 

, such a spacetime acquires a constant asymptotic speed at large distances. This speed, 
calibrated empirically, is proportional to the fourth root of the mass of the central body 
in compliance with the Tully-Fisher relation. A variance of MOND's gravity emerges as 
\ an inevitable consequence of the proposed formalism. It has also been shown (Mendoza 
^p" 1 . et al. 2006) that a) the gravitational waves in this spacetime propagate with the speed 
2 ■ of light in vacuum and b) there is a lensing effect added to what one finds in the classic 
GR. 

03 



1. Introduction 

Convinced of cosmic speed up and not finding the dark energy hypotheses a compelling 
explanation, some cosmologists have looked for alternatives to Einstein's gravitation (Def- 
fayet et al. 2002; Freese et al. 2002; Ahmed et al. 2002; Dvali et al. 2003; Capozziello et al. 
2003; Carroll et al. 2003; Norjiri et al. 2003, 2004, and 2006; Das et al. 2005; Sotiriou 2005; 
Woodard 2006). There is a parallel situation in galactic studies. Dark matter hypotheses, 
intended to explain the flat rotation curves of spirals or the large velocity dispersions in 
ellipticals, have raised more questions than answers. 

Alternatives to newtonian dynamics have been proposed but have had their own 
critics. Foremost among such theories, the Modified Newtonian Dynamics (MOND) of 
Milgrom (1983 a,b,c) is able to explain the flat rotation curves (Sandres et al. 1998 
and 2002) and justify the Tully-Fisher relation with considerable success. But it is often 
criticized for the lack of an axiomatic foundation; see, however, Bekenstein's (2004) TeVeS 
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theory where he attempts to provide such a foundation by introducing a tensor, a vector, 
and a scalar field into the field equations of GR. 

Here we are concerned with galactic problems. We suggest following cosmologists and 
look for a modified Einstein gravity tailored to galactic environments. In Sects. 2 and 3 
we design an action integral, different but close to that of Einstein-Hilbert, and find a 
spherically symmetric static solution to it. In Sect. 4 we analyze the orbits of test objects 
moving in this modified spacetime and demonstrate the kinship of the obtained dynamics 
with MOND. Section 5 is devoted to concluding remarks. 



2. A modified field equation 



The model we consider is an isolated mass point. As an alternative to the Einstein-Hilbert 
action, we assume 
1 



■ s = 2 / f(R)V=9d 4 x, 



(1) 



where R is the Ricci scalar and f(R) an as yet unspecified, but differentiable function of 
R. Variations in S with respect to the metric tensor lead to the following field equation 
(Capozziello et al. 2003): 



r I t - (h -h A x 1 



h' 
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where h = df /dR. The case f(R) = R + constant and h = 1 gives the Einstein field 
equation with a cosmological constant included in it. For the purpose of galactic studies, 
we envisage a spherically symmetric static Schwarzschild-like metric, 



ds 2 = -B(r)dt 2 + A{r)dr 2 + r 2 ( 
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Equation ([J| is the combination R tt /B + R rr /A, Eq. ((SJ) is R rr /A — Ree/r 2 , and Eq. © 
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is the rr-component of the field equation. Finally, Eq. ((7]) is from the contraction of Eq. 
(J5J. In principle, for a given h (or/) one should be able to solve the four Eqs. I©-© for 
the four unknowns, A, B, R, and / (or h), as functions of r. 



3. Solutions of Equations (4)-(7) 



We are interested in solutions that differ from those of the classic GR by small amounts. 
For the classic GR one has h = 1 and A{r)B(r) = 1. Here, we argue that, if the combi- 
nation B'/B + A'/A is a well-behaved differential expression, it should have a solution of 
the form A(r)B{r) — g(r). Furthermore, g(r) should differ from 1 only slightly, in order 
to remain in the vicinity of GR. There are a host of possibilities. For the sake of argument 
let us assume g(r) = (r/s) a fall cvln(r/s), where a is a small dimensionless parameter 
and s is a length scale of the system to be identified shortly. Equation (J4J) splits into 



B A a .„ /r\ Q 

h -r = -, AB = - 

B A r \sJ 

h" - -~h! + ^h = 0. 



(8) 
(9) 

^ Of these, 



2 r 

Equation ([9|) has the solution h = (r/s)' 3 , — a + O (a 2 ), and 1 — ^a + O [a 1 - 
the solution h « {r/s) a satisfies the requirement h — > 1 as a — > 0. The second solution 
is discarded. Substituting AB = h = (r/s) a in Eq. (O gives 



1 



1 

A ~ (1 - a) 
r\ a 1 
s 



/s\( 1 - a / 2 ) , /r\ 



f,2(l-a/2) 



A 



(10) 
(11) 



where A is a constant of integration. Actually there is another constant of integration 
multiplying the (s/r) term. We have, however, absorbed it in the expression for s that we 
now define. For a = 0, Eqs. pTJ|) and (fTTj) are recognized as the Schwarzschild-deSitter 
metric. Therefore, s is identified with the Schwarzschild radius of a central body, 2GM/c 2 , 
and A with a dimensionless cosmological constant. Substitution of Eqs. (|10p and pip into 
Eqs. © and © gives 



/ = 



R = 



1 



(1 — a) r 2 

3 1 

(1 — a) r 2 



+ (4-a)A(£) 
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The Ricci scalar of the Schwarzschild space is zero and that of the deSitter or the 
Schwarzschild-deSitter space is constant. For non zero cv, however, R is somewhere be- 
tween these two extremes. At small distances it increases as r~ 2 and at large r's it behaves 
as s~ 2 (s/r) a sw s~ 2 (l — alnr/s). The spacetime is asymptotically neither flat nor deSit- 
terian. Cosmologists may find this variable Ricci scalar relevant to their purpose ( see 
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also Brevik et al, 2004, for a different modification of the Schwarzschild-deSitter metric). 
Likewise, we began with / as a function of R rather than r. Elimination of r between 
Eqs. (fl2|) and (|13[) provides one in terms of the other. For A = 0, one easily finds 

/ = (3a) a/2 s- a iJ^°/ 2 ) w R[l - | ln(s 2 i?) + | ln(3a)]. (14) 

Once more we observe the mild logarithmic correction to the classic GR. 

4. Applications to galactic environments 

In this section we demonstrate that 

— The logarithmic modification of the Einstein-Hilbert action, in the weak field regime, 
results in a logarithmic correction to the newtonian potential. A test star moving in 
such a potential acquires a constant asymptotic speed, Woo = cy/a/2. 

— The asymptotic speed cannot be independent of the central mass. We resort to the 
observed rotation curves of spirals to find this dependence. 

— The high- and low- acceleration limits of the weak-field regime are the same as those 
of MOND. A kinship with MOND follows. 

4-1- Orbits in the spacetime of Equations 

We assume a test star orbiting a central body specified by its Schwarzschild radius, 
2GM/c 2 . We choose the orbit in the plane 9 = ir/2. The geodesic equations for r, p, and 
t are 

d 2 r 1A' fdrV r fdip\ 2 1 B> ( dt\ 2 , . 

d^ + 2A{Tr) -Ail) + 2tU) =°' (15) 

dip\ 1 d 2 (p 2 dr 



dr J dr 2 r dr ' ^ ^ 



dt V 1 d 2 t B' dr 
Jt) dT 2+ ~BdT' 
respectively. Equations p^|) and (jTTJ) immediately integrate into 



) — + — — = ' ( 17 ) 



r dip/dr — J, a constant, (18) 



dt/dr = l/B. (19) 

Substituting the latter into Eq. (|15[) and assuming a circular orbit, dr/dr = 0, gives 
J 2 _ 1 B' _ 1 /r\°> B' 

^^2B^^2\l) S 1 ' (20) 
where we have used Eq. Ijlip to eliminate A. In galactic environments what one measures 
as the circular orbital speed is 

rdip r dip dr \J~B J 



v^Bdt VB dr dt 



(21) 
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Eliminating J between Eqs. ([2Tjl and (|20|) gives 

(22) 



2 1 rB ' 1 

B = 2~2 



Further substitution for B from Eqs. (fTTj) and (fTU|) yields 

,(!-"/2) , o\ /r\ 2 ( 1 -"/2)l 



1 1 



(£)^-"+2A(5) 



w = o« + o (1 - «/2) r L ,., (1 - a/2) Wr , 2( i- a/2 )T ( 23 ) 



2 2 



To put Eq. ([23]) in a tractable form: 

— We neglect the A term and substitute s = 2GM/c 2 . 

— We adopt the approximation x~ a — exp(— a\nx) = 1— alnx + O (a 2 ). 

— The terms containing s are small. We retain only the first order terms in s. 

— v is measured in units of c. We restore it hereafter. 



With these provisions, Eq. (|23p reduces to 

1 ' (24) 



2 1 iGM 

v = —ac H 

2 r 



1 f / 2GM \ 
1 - -a { 1 + In 
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A plot of v 2 as a function of r has the horizontal asymptote \ol<? 



4-2. Determination of a 



The asymptote in Eq. (|24|) cannot be a universal constant. It is not possible to imagine 
that a galaxy and a speck of dust dictate the same speed for distant passing objects. The 
parameter a should depend on the mass of the gravitating body residing at the origin, 
because any localized matter will betray no characteristics other than its mass when 
sensed from far distances. To find the mass dependence of a we resort to observations. 
From Sanders and Verheijn (1998) and Sanders and McGaugh (2002), we have compiled 
a list of 31 spirals for which total masses, asymptotic orbital speeds, and velocity curves 
are reported. The figures in their papers contain the observed circular speeds and the 
newtonian ones derived from the observed mass of the stellar and HI components of the 
galaxies. We have selected those objects that a) have a noticeable horizontal asymptote, 
b) have fairly reduced newtonian speeds by the time the flat asymptote is approached, 
and c) do not possess anomalously high HI content to hinder estimates of the total mass 
and the size of the galaxy. We also made the assumption that the total HI and stellar 
mass are distributed spherically symmetrically and mimic a point mass if observed from 
far distances. The relevant data along with a = 2v 2 a /c 2 are reported in the table, and 
the figure is a log-log plot of the calculated a versus the mass. A power law fit to the 
data gives 

a = (3.07 ± 0.18) x 10- 7 (Af/10 10 M o ) a494 . (25) 
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It is important to note that Eq. (f2"5"l) is not a consequence of the present theory, but 
rather an empirical relation dictated by observations and based on the masses and the 
asymptotic speeds of a selected list of galaxies reported by Sanders et al. Together with 
the popularly accepted rule that the masses and the luminosities of spirals are linearly 
related, it leads to a Tully-Fisher (TF) relation, Luminosity oc Woo 4 ' 05 - Observational 
actualities, however, are complicated. In a recent paper, Kregel et al. (2005) distinguish 
between different TF relations based on the luminosity, disk mass, maximum disk stel- 
lar mass, baryonic mass (meaning stellar+HI mass), baryonic + bulge mass, etc. The 
reported exponents range from 3.23 ± 0.36 to 4.2 ± 0.23, depending on the type of qual- 
ification; see also Gurovich et al. (2004). A more elaborate discussion of the issue falls 
beyond the scope of the present paper. 

The main sources of error in Eq. (|25|) . both in the exponent and in the slope, are 
a) the estimates of the total masses of the galaxies, b) the judgment whether what one 
measures as the asymptotic speed is indeed the orbital speed at the far outskirts of the 
galaxy, c) the popular assumption that the masses and luminosities of the spirals are 
linearly related, and finally, d) our heuristic assumption that the galaxies can be treated 
as spherically symmetric objects. In spite of all these uncertainties, we note that the 
exponent 0.494 is astonishingly close to 0.5, the figure that one finds from MOND. We 
also demonstrate in the following section that the slope 3.00 x 10~ 7 is in very good 
agreement with the characteristic acceleration of MOND. 




-1 -0.5 0.5 1 1.5 

x=\oa (M/10 10 M ) 

Fig. 1. A log- log plot of a versus M. The equation for the power law fit is shown in the 
legend. 

4.3. Kinship with MOND 

We recall that in the weak-field approximation, the newtonian dynamics is derived from 
the Einsteinian one by writing the metric coefficient B = (l + 20/c 2 ) , (j> = GM/r and by 
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Table 1. The data in the first four columns are from Sanders et al. 2002. The last two 
columns show the empirical relation between the asymptotic speeds and the masses of 
the galaxies. 



Galaxy 


H 


M 


Voo 


2{Voo/C) 


ao 




kpc 


1U Mq 


km/s 


X 1U 


X 1U 


NGC 5533 


72.0 


22.0 


250 


13.9 


3.02 


NGC 3992 


30.0 


16.22 


242 


13.0 


3.28 


NGC 5907 


32.0 


10.8 


214 


10.2 


3.15 


NGC 2998 


48.0 


11.3 


213 


10.1 


3.05 


NGC 801 


60.0 


12.9 


218 


10.6 


3.00 


NGC 5371 


40.0 


12.5 


208 


9.61 


2.76 


NGC 4157 


26.0 


5.62 


185 


7.61 


3.24 


NGC 4217 


14.5 


4.50 


178 


7.04 


3.35 


NGC 4013 


27.0 


4.84 


177 


6.96 


3.19 


NGC 4088 


18.8 


4.09 


173 


6.65 


3.32 


NGC 4100 


19.8 


4.62 


164 


5.98 


2.81 


NGC 3726 


28.0 


3.24 


162 


5.83 


3.26 


NGC 4051 


10.6 


3.29 


159 


5.62 


3.12 


NGC 4138 


13.0 


3.01 


147 


4.82 


2.80 


NGC 2403 


19.0 


1.57 


134 


3.99 


3.19 


UGC 128 


40.0 


1.48 


131 


3.81 


3.14 


NGC 3769 


33.0 


1.33 


122 


3.31 


2.88 


NGC 6503 


21.8 


1.07 


121 


3.25 


3.14 


NGC 4183 


18.0 


0.93 


112 


2.79 


2.89 


UGC 6917 


9.0 


0.74 


110 


2.69 


3.12 


UGC 6930 


14.5 


0.73 


110 


2.69 


3.14 


M 33 


9.0 


0.61 


107 


2.54 


3.24 


UGC 6983 


13.8 


0.86 


107 


2.54 


2.74 


NGC 7793 


6.8 


0.51 


100 


2.22 


3.10 


NGC 300 


12.4 


0.35 


90 


1.80 


3.02 


NGC 5585 


12.0 


0.37 


90 


1.80 


2.94 


NGC 6399 


6.8 


0.28 


88 


1.72 


3.22 


NGC 55 


10.0 


0.23 


86 


1.64 


3.39 


UGC 6667 


6.8 


0.33 


86 


1.64 


2.83 


UGC 6923 


4.5 


0.24 


81 


1.46 


2.95 


UGC 6818 


6.0 


0.14 


7o 


1.18 


3.12 



R: radius of the galaxy (kpc); M: stellar + HI mass (10 Mq); u^: asymptotic speed (km/sec); 
a : 2(v oo /c) 2 M- 0A94 . 
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expanding all relevant functions and equations up to the first order in </>/c 2 . fn a similar 
way one may find our modified newtonian dynamics from the presently modified GR by 
expanding B(r) of Eq. (|Iip up the first order in a and s/r. Thus 

B(r) = 1 + a + ahx(r/s) - s/r = I + 2</>(r)/c 2 , (26) 

where the second equality defines 4>{r). Let us write Eq. (j25|) (with slight tolerance) as 
a = «o(GM / 'GMq) 1 / 2 and find the gravitational acceleration 



g = \d<f>/dr\ = (a Q g n ) 1/2 +g n (27) 

= g n for g n > a Q 

= (a g n ) 1/2 for a > g n -> 0, 

where we have denoted 

a = alc A /AGM Q and g n = GM/r 2 . (28) 

The limiting behaviors of g are the same as those of MOND. One may, therefore, com- 
fortably identify ao as MOND's characteristic acceleration and calculate cto anew from 
Eq d25]). For a = 1.2 x f(T 8 cm/sec 2 , one finds 

a = 2.8 x f(T 12 (M/M ) 1/2 . (29) 

It is gratifying how close this value of a is to the one in Eq. (|25p and how similar 
the low and high acceleration limits of MOND and the present formalism are, in spite of 
their totally different and independent starting points. It should also be noted that there 
is no counterpart to the interpolating function of MOND here. 

5. Concluding remarks 

We have developed an f(R) oc R}~ a ' 2 gravitation that is essentially a logarithmic mod- 
ification of the Einstein-Hilbert action. In spherically-symmetric static situations, the 
theory allows a modified Schwarzschild-deSitter metric. This metric in the limit of weak 
fields gives a logarithmic correction to the newtonian potential. From the observed asymp- 
totic speeds and masses of spirals we learn that the correction is proportional to almost 
the square root of the mass of the central body. Flat rotation curves, the Tully-Fishcr 
relation (admittedly with some reservations) , and a version of MOND emerge as natural 
consequences of the theory. 

Actions are ordinarily form invariant under the changes in sources. Mass dependence 
of a destroys this feature and any claim for the action-based theory should be qualified 
with such reservation in mind. This, however, should not be surprising, for it is understood 
that all alternative gravitations, one way or another, go beyond the classic GR. One 
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should not be surprised if some of the commonly accepted notions require re-thinking 
and generalizations. 

Since the appearance of an earlier version of this paper in arXiv, Mendoza et al.(2006) 
have investigated the gravitational waves and lensing effects in the proposed spacetime. 
They find the following: a) in any f(R) = R n gravitation, gravitational waves travel with 
the speed of light in a vacuum, and b) in the present spacetime, there is a lensing in 
addition to what one finds in the classical GR. Their ratio of the additional deflection 
angle of a light ray, 5(3, to that in GR, Pgr, can be reduced to 



where r m is the impact parameter of the impinging light. The proportionality of 5(3 to 
a is expected, because the proposed metric is in the neighborhood of GR. Its increase 
with increasing r m also should not be surprising, since the theory is designed to highlight 
unexpected features at far rather than nearby distances. 

Soussa et al. (2004) maintain that "no purely metric-based, relativistic formulation 
of MOND whose energy functional is stable can be consistent with the observed amount 
of gravitational lensing from galaxies" . For at least two reasons, this no-go theorem does 
not apply to what we have highlighted above as the kinship with MOND: 

a) Apart from their common low and high acceleration regimes, the two theories are 
fundamentally different. The gravitational acceleration in the weak field limit of the 
present theory is the newtonian one to which a small \jr correction is added. That of 
MOND, on the other hand, is a highly nonlinear function of the newtonian acceleration 
through an arbitrary interpolating function. 

b) More important, however, is one of the authors' assumptions that "the gravitational 
force is carried by the metric, and its source is the usual stress tensor" . This is not the 
case in the present theory. Although we have only worked out the vacuum solution for 
a point source, the mass dependence of the exponent a in Eqs. (|10p and (|lip makes the 
theory different from what the assumption requires. 

There are two practices for obtaining the field equations of /(i?) gravity, the metric 
approach, where g^s are considered as dynamical variables, and that of Palatini, where 
the metric and the affine connections are treated as such (see Magnano 1995 for a review). 
Unless f(R) is linear in R, the resulting field equations are not identical (see Ferraris et 
al. 1994). The metric approach is often avoided for leading to fourth-order differential 
equations. It is also believed to have instabilities in the weak field approximations (see 
e. g., Sotiriou 2005 and also Amarzguioui et al. 2005). In the present paper we do not 
initially specify f(R)- Instead, at some intermediate stage in the analysis we adopt an 
ansatz for df(R)/dR as a function of r and work from there to obtain the metric, the 



5(3/(3 G r 




(30) 
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Ricci scalar, and eventually f(R). This enables us to avoid the fourth-order equations. 
This trick should work in other contexts, such as cosmological ones. 

The theory presented here is preliminary. Further investigations are needed from 
both formal and astrophysical points of view. The author's list of priorities include the 
following: 

— Stability of the metric of Eqs. (JTUJ) fc fTTj) . The approach should be to impose a small 
perturbation 5g^ v on the metric, linearize the field Eq. ([2]), and ask for the condition 
of stability of the metric. Such a condition, if it exists at all, might throw some 
light on the mass dependence of a, the empirical relation of Eq. (|29|) . Managing the 
linear problems is straightforward. Here, however, the bookkeeping is extensive and 
laborious. 

— Extension of the theory, at least in the weak field regime, to many body systems and 
to cases with a continuous distribution of matter, in order to obtain the metric inside 
the matter. 

— Developing the theory beyond the first order in a 

— Solar system tests of the theory. 

— Possible cosmological implications of the theory. 

Acknowledgement: The author wishes to thank Bahram Mashhoun and Naresh Dad- 
hich for comments and helpful suggestions. Reza Saffari has pointed out a typographical 
error in the earlier version of the paper, Eq. (15). The error had propagated making the 
sign of a appear positive in the factor multiplying the term GM/r in Eq. (24). This is 
corrected here. 
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